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Abstract. Effects of lattice-QCD-inspired equations of state and continuous emission on some 
observables are discussed, by solving a 3D hydrodynamics. The particle multiplicity as well as 
V2 are found to increase in the mid-rapidity. We also discuss the effects of the initial-condition 
fluctuations. 

Keywords: LatticeQCD equations of state, hydrodynamic model 
PACS: 24.10.Nz,25.75.-q,25.75.Ld 

1. HYDRODYNAMIC MODELS 

Hydrodynamics is one of the main tools for studying the collective flow in high-energy 
nuclear collisions. Here, we shall examine some of the main ingredients of such a de- 
scription and see how likely more realistic treatment of these elements may affect some 
of the observable quantities. The main components of any hydrodynamic model are the 
initial conditions, the equations of motion, equations of state and some decoupling pre- 
scription. We shall discuss how these elements are chosen in our studies. 

Initial Conditions: In usual hydrodynamic approach, one assumes some highly sym- 
metric and smooth initial conditions (IC). However, since our systems are small, large 
event-by-event fluctuations are expected in real collisions, so this effect should be taken 
into account. We introduce such IC fluctuations by using an event simulator. As an ex- 
ample, we show here the energy density for central Au+Au collisions at 130A GeV, 




FIGURE 1. The initial energy density at 77 = is plotted in units of GeV/fm 3 . One random event is 
shown vs. average over 30 random events (~ smooth initial conditions in the usual hydro approach). 



given by NeXuS 1 01- Some consequences of such fluctuations have been discussed 
elsewhere|H,0,|7|]. We shall discuss some others in Sec|^l 

Equations of Motion: In hydrodynamics, the flow is governed by the continuity equa- 
tions expressing the conservation of energy-momentum, baryon-number and other con- 
served charges. Here, for simplicity, we shall consider only the energy-momentum and 
the baryon number. Since our systems have no symmetry as discussed above, we de- 
veloped a special numerical code called SPheRIO (Smoothed Particle hydrodynamic 
evolution of Relativistic heavy IOn collisions) based on the so called Smoothed- 
Paricle Hydrodynamics (SPH) algorithm [9]. The main characteristic of SPH is the 
parametrization of the flow in terms of discrete Lagrangian coordinates attached to small 
volumes (called "particles") with some conserved quantities. 

Equations of State: In high-energy collisions, one often uses equations of state (EoS) 
with a first-order phase transition, connecting a high-temperature QGP phase with a low- 
temperature hadron phase. A detailed account of such EoS may be found, for instance, in 
0]. We shall denote them 10PT EoS. However, lattice QCD showed that the transition 
line has a critical end point and for small net baryon surplus the transition is of crossover 
type ifioll . The following parametrization may reproduce this behavior, in practice: 

P = XP h + (1-X)P q + 25/^(Pq-P h ) 2 + 48 , (1) 
s = Xs h + {\-X)sqi (2) 

e = Xe H + (l-X)£ Q -2 [l + (m/m,) 2 ] S/ ^(P Q -P H ) 2 + ^8 , (3) 
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FIGURE 2. A comparison of e(T), s(T) and P(T) as given by our parametrization with a critical point 
(solid lines) and those with a first-order phase transition (dashed lines). 



1 Many other simulators, based on microscopic models, e.g. HIJING \2], VNI H, URASiMA 14], • • ■, 
show such event-by-event fluctuations. 




FIGURE 3. Plots of s/e and P as function of e for the two EoS shown in Figure|2] 



where A = [1 - (P Q - Ph)/ \J(Pq-Ph) 2 + ]/2 and suffixes Q and H denote those 
quantities given by the MIT bag model and the hadronic resonance gas, respectively, 
and 8 = 5(jU&) = <5oexp[— (jU^/jU c ) 2 ], with \i c =const. As is seen, when 8(flb) ^ 0, the 
transition from hadron phase to QGP is smooth. We could choose 8 (fib) so to make it 
exactly when fib > jU c , to guarantee the first-order phase transition there. However, in 
practice our choice above showed to be enough. We shall denote the EoS given above, 
with 8q 7^ 0, CP EoS. Let us compare, in Figure |2j e(T), s(T) and P(T), given by the 
two sets of EoS. one can see that the crossover behavior is correctly reproduced by our 
parametrization for CP EoS, while finite jumps in e and s are exhibited by 10PT EoS, 
at the transition temperature. It is also seen, as mentioned above, that at fib ~ 0.4 GeV 
the two EoS are indistinguishable. Now, since in a real collision what is directly given is 
the energy distribution at a certain initial time (besides n^, s, etc.), whereas T is defined 
with the use of the former, we plotted some quantities as function of e in Figure |3] 
One immediately sees there some remarkable differences between the two sets of EoS: 
naturally p is not constant for CP EoS in the crossover region; moreover, s is larger. We 
will see in Sec|2]fhat these features affect the observables in non-negligible way. 

Decoupling Prescription: Usually, one assumes decoupling on a sharply defined hy- 
persurface. We call this Sudden Freeze Out (FO). However, since our systems are small, 
particles may escape from a layer with thickness comparable with the systems' sizes. 
We proposed an alternative description called Continuous Emission (CE) 111 ill which, 
as compared to FO, we believe closer to what happens in the actual collisions. In CE, 
particles escape from any space-time point jc^, according to a momentum-dependent 
escaping probability &(x,k) = exp [— p(x') ov dV] . To implement CE in SPheRIO 
code, we had to approximate it to make the computation practicable. We took & on the 
average, i.e., 

&>(x,k) -> (^(x,k)} = &>{x) = exp (-K s 2 /\ds/dx\) . (4) 
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FIGURE 4. 77 and p T distributions for the most central Au+Au at 200A GeV. Results of CP EoS and 
10PT EoS are compared. The data are from PHOBOS Collab.fll. 



The last equality has been obtained by making a linear approximation of the density 
p(x') = ocs(x') and k = 0.5 a (ov) is estimated to be 0.3, corresponding to (av) ~ 
2 fm 2 . It will be shown in Sec.|2]that CE gives important changes in some observables. 



2. RESULTS 

Let us now show results of computation of some observables, as described above, for 
Au+Au at 200A GeV. We start computing r\ and pj distributions for charged particles, 
to fix the parameters. Then, V2 and HBT radii are computed free of parameters. 

Pseudo-rapidity distribution: Figure |3] shows that the inclusion of a critical end point 
increases the entropy per energy. This means that, given the same total energy, CP 
EoS produces larger multiplicity, which is clearly shown in the left panel of Figure |4j 
especially in the mid-rapidity region. Now, we shall mention that, once the equations of 
state are chosen, fluctuating IC produce smaller multiplicity, for the same decoupling 
prescription, as compared with the case of smooth averaged IC CD- 




FIGURE 5. Left: tj distribution of V2 for charged particles in the centrality (15-25)% Au+Au at 200A 
GeV, computed with fluctuating IC. The vertical bars indicate dispersions. The data are from PHOBOS 
Collab.| 13]. Right: v 2 distribution in the interval 0.48 < 77 < 0.95, corresponding to CP EoS and CE. 
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FIGURE 6. kj dependence of HBT radius Rl for n in the most central Au+Au at 200A GeV, computed 
with fluctuating IC. The data are from PHENIX Collab. fil . 



Transverse-Momentum Distribution: As discussed in Sec. [H since the pressure does 
not remain constant in the crossover region, we expect that the transverse acceleration 
is larger for CP EoS, as compared with 10PT EoS case. In effect, the right panel of 
Figure |4| does show that pj distribution is flatter for CP EoS, but the difference is 
small. The freezeout temperature suggested by r\ and pj distributions turned out to 
beT f ~ 135-140MeV. 

Elliptic-Flow Parameter V2'. We show, in Figure El results for the r\ distribution of V2 
for Au+Au collisions at 200A GeV. As seen, CP EoS gives larger V2 , as a consequence 
of larger acceleration in this case as discussed in Sec[TJ Notice that CE makes the curves 
narrower, as a consequence of earlier emission of particles, so with smaller acceleration, 
at large- 1 T] | regions. Due to the IC fluctuations, the resulting fluctuations of V2 are large, 
as seen in Figures |5J It would be nice to measure such a V2 distribution, which would 
discriminate among several microscopic models for the initial stage of nuclear collisions. 

HBT Radii: Here, we show our results for the HBT radii, in Gaussian approximation as 
often used, for the most central Au+Au collisions at 200A GeV. As seen in Figures |()] and 




FIGURE 7. kj dependence of HBT radii R s and R for pions in the most central Au+Au at 200A GeV, 
computed with event-by-event fluctuating IC. The data are from PHENIX Collab.| 14]. 



|7J the differences between CP EoS results and those for 10PT EoS are small. For R s , 
and especially for R D , one sees that CP EoS combined with continuous emission gives 
steeper kj dependence, closer to the data. However, there is still numerical discrepancy 
in this case. 

3. CONCLUSIONS AND OUTLOOKS 

In this work, we introduced a parametrization of lattice-QCD EoS, with a first-order 
phase transition at large jU& and a crossover behavior at smaller \X\, . By solving the hy- 
drodynamic equations, we studied the effects of such EoS and the continuous emission. 
Some conclusions are: i) The multiplicity increases for these EoS in the mid-rapidity; ii) 
The px distribution becomes flatter, although the difference is small; Hi) V2 increases; 
CE makes the r\ distribution narrower; iv) HBT radii slightly closer to data. 

In our calculations, the effect of the continuous emission on the interacting component 
has not been taken into account. A more realistic treatment of this effect probably makes 
R smaller, since the duration for particle emission becomes smaller in this case. Another 
improvement we should make is the approximations we used for &>(x,p). 
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